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Abstract 

The quantum numbers J^*^ = 1^^ of the X{3872) and the proximity of its mass to the D*^D^ 
threshold imply that it is either a loosely-bound hadronic molecule whose constituents are a su- 
perposition of D*°D° and L»°Z)*° or it is a virtual state of charm mesons. The line shapes of the 
X(3872) can discriminate between these two possibilities. At energies within a few MeV of the 
]j*0]jO threshold, the lines shapes of the X produced in B ^ K transitions are determined by its 
binding energy and its width. Their normalizations are determined by a short-distance constant 
that is different for and B^ — > K^. At energies comparable to the 8 MeV splitting be- 

tween the D*^D^ and D*^D~ thresholds, the charged meson channels D*~^D~ and D'^D*~ have 
a significant effect on the line shapes of the X. We calculate the line shapes taking into account 
the resonant coupling between the charged and neutral 1"^^ channels. The line shapes and their 
normalizations depend on one additional scattering parameter and two additional short-distance 
constants associated with the B ^ K transitions. The line shapes of the X resonance depend on 
its decay channel; they are different for J / ijj tt^ -k^ , J / ip ir'^ Tr~ ir^ , and D^D^tt^. The line shapes are 
also different for X produced in B~^ decays and in B^ decays. Some conceptual errors in previous 
work on this problem are pointed out. 

PACS numbers: 12.38.-t, 12.39. St, 13.20.Gd, 14.40. Gx 
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I. INTRODUCTION 



The X(3872) is a cc resonance near 3872 MeV discovered in 2003 by the Belle Collabo- 
ration [1] and subsequently observed by the CDF, Babar, and DO Collaborations [2-4]. In 
addition to the discovery decay mode ,//?/-' tt+tt", the X has been observed to decay into 
J/'ip^, J / ip TT~^ IT' TT^ , [5] and D^D^ir^ [6, 7]. The decay into J/ip'~f implies that the X is even 
under charge conjugation. An analysis by the Belle Collaboration of the decays of X into 
J/ipTi~^Ti~ strongly favors the quantum numbers J^*^ = l"*""*", but does not exclude 2"'""'" [8]. 
An analysis by the CDF Collaboration of the decays of X into J/t/jn'^n' is compatible with 
the Belle constraints [9]. The tiny phase space available for the decay into D^D^-k^ rules 
out J = 2, leaving 1++ as the only option. 

An important feature of the X(3872) is that its mass Mx is extremely close to the 
]j*o£)0 threshold. The PDG value for Mx comes from combining measurements of X in 
the J/il)TT^TT~ decay mode [10]. After taking into account a recent precision measurement 
of the mass by the CLEO Collaboration [11], the difference between the PDG value for 
Mx and the D*^D^ threshold is 

Mx - (M,o + Mo) = -0.6 ± 0.6 MeV, (1) 

where M^o and Mq are the masses of D*^ and . The negative central value in Eq. (1) is 
compatible with the X being a bound state of the charm mesons. The measured mass of the 
near-threshold enhancement in D°5°7r° is about 4 MeV above the D*^D^ threshold [6, 7]. 
This value is compatible with X being a virtual state of charm mesons. It differs from the 
mass in Eq. (1) by more than two standard deviations, which raises the question of whether 
the decays into J/ipTi^n' and D^D^ir^ are coming from the same resonance. 

The proximity of the mass of the X(3872) to the D*^D^ threshold has motivated its 
identification as a weakly-bound molecule whose constituents are a superposition of the 
charm meson pairs D*^D^ and D^D*^ [12-15]. The establishment of the quantum numbers 
of the X(3872) as 1++ makes this conclusion almost unavoidable. The reason is that these 
quantum numbers allow S-wave couplings of the X to D*^D^ and D^D*^. Nonrelativistic 
quantum mechanics implies that a resonance in an S-wave channel near a 2-particle threshold 
has special universal features [16]. Because of the small energy gap between the resonance 
and the 2-particle threshold, there is a strong coupling between the resonance and the 
two particles. This strong coupling generates dynamically a large length scale that can be 
identified with the absolute value of the S-wave scattering length a of the two particles. 
Independent of the original mechanism for the resonance, the strong coupling transforms 
the resonance into a bound state just below the two-particle threshold if a > or into a 
virtual state just above the two-particle threshold if a < 0. If a > 0, the bound state has a 
molecular structure, with the particles having a large mean separation of order a. 

To see that the universal features of an S-wave threshold resonance are relevant to the 
X(3872), we need only note that its binding energy is small compared to the natural energy 
scale associated with pion exchange [17]: m^/(2M*oo) ~ 10 MeV, where M*oo is the reduced 
mass of the two constituents. The universal features of the X(3872) were first exploited by 
Voloshin to describe its decays into D°£)°7r° and D^D^^, which can proceed through decay 
of the constituent D*^ or D*^ [15]. Universahty has also been apphed to the production 
process B KX [18, 19], to the line shapes of the X [20], and to decays of X into J/il) 
and pions [21]. These applications rely on factorization formulas that separate the length 
scale a from all the shorter distance scales of QCD [20] . The factorization formulas can be 
derived using the operator product expansion for a low-energy effective field theory [22] . 
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Other interpretations of the X(3872) besides a charm meson molecule or a charm meson 
virtual state have been proposed, including a P-wave charmonium state or a tetraquark 
state. (For a review, see Ref. [23].) If the charmonium or tetraquark models were extended 
to include the couphng of the X to D*^D^ and D^D*^ scattering states, the universal features 
of an S-wave threshold resonance imply that the tuning of the binding energy to the threshold 
region would transform the state into a charm meson molecule or a virtual state of charm 
mesons. Any model of the X(3872) that does not take into account its strong coupling to 
charm meson scattering states should not be taken seriously. 

Given that the quantum numbers of the X{3S72) are l"'""'", the measured mass Mx in 
Eq. (1) implies unambiguously that X must be either a charm meson molecule or a virtual 
state of charm mesons. The remaining challenge is to discriminate between these two pos- 
sibilities. If the X was sufficiently narrow, there would be clear qualitative differences in 
its fine shapes between these two possibilities. We first consider the D^D'^n^ decay mode, 
which has a contribution from the decay of a constituent If the X was a charm meson 
molecule, its line shape in D^D^n^ would consist of a Breit-Wigner resonance below the 
]j*0£)0 threshold and a threshold enhancement above the D*^D^ threshold. If the X was a 
virtual state, there would only be the threshold enhancement above the D*^D^ threshold. 
We next consider decay modes that have no contributions from the decay of a constituent 
D*^, such as J / ip Ti^ Ti~ . If the X was a charm meson molecule, its line shape in such a decay 
mode would be a Breit-Wigner resonance below the threshold. If the X was a virtual 

state, there would only be a cusp at the D*^D^ threshold. The possibility of interpreting 
the X(3872) as a cusp at the D*^D^ threshold has been suggested by Bugg [24]. Increasing 
the width of the X provides additional smearing of the line shapes. This makes the quali- 
tative difference between the line shapes of a charm meson molecule and a virtual state less 
dramatic. To discriminate between these two possibilities therefore requires a quantitative 
analysis. 

There have been two recent analyses of data on B ^ K -\- X that shed fight on the issue 

of whether the X is a bound state or a virtual state. Hanhart et al. [25] analyzed the data 
on B+ ^ K+ + J I'll) TT+TJ- and B+ ^ K+ + D^D^n^ from the Belle and Babar Collabora- 
tions using a model for the scattering amplitude in the D*^D^'' + D'^D*'^ channel that is a 
generafization of the Flatte parametrization for a near-threshold resonance. They concluded 
that the D^D^n^ threshold enhancement observed by the Belle and Babar Collaborations 
is compatible with the X(3872) only if the X is a virtual state. One flaw in the analysis of 
Ref. [25] is that it did not take into account the width of the constituent D*^. They also 
assumed incorrectly that a bound state below the D*°D^ threshold would not decay into 

In Ref. [26], we derived the line shapes of the X(3872) near the D*^D^ threshold from 
the assumption of an S-wave resonance in the neutral charm meson channel D*^ D'^ + D*^ . 
We developed expressions for the line shapes that take into account the width of the 
meson and the inelastic scattering channels of the charm mesons. An analysis of the data on 
B+ ^ K+ + J/tlj TT+vr- and B+ ^ K+ + D^D^tc^ from the Belle Collaboration indicated that 
the data preferred the X(3872) to be a bound state but a virtual state was not excluded. The 
most important lesson of the analyses of Refs. [25] and [26] is that the measured difference 
between the masses of the X(3872) in the J/V' tt+tt" and D^D^n'^ decay channels is consistent 
with it being a charm meson molecule or a virtual state of charm mesons. 

In this paper, we generalize the results of Ref. [26] for the line shapes of the X(3872) to 
take into account the resonant coupling between the neutral charm meson channel and the 
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charged charm meson channel D*^D^ + D~^D*~ . In Sec. II, we summarize the results of 
Ref. [26] for the energy-dependent widths of virtual D* mesons. In Ref. [26], we developed an 
expression for the resonant scattering amplitude for the neutral charm meson channel that 
takes into account the D*° width and inelastic charm meson scattering channels. In Sec. Ill, 
we extend that result to the three scattering amplitudes for the resonantly coupled neutral 
and charged charm meson channels. In Ref. [26], we derived factorization formulas for the 
line shapes of X(3872) in the decays B ^ K + X that take into account the resonance in the 
neutral charm meson channel. In Section IV, we extend those results to take into account 
the resonant coupling to the charged charm meson channel. In Sec. V, we summarize our 
results. 



II. MASSES AND D* WIDTHS 

When we consider the decays of the D* mesons, there are particles with six different 
masses that must be considered. We therefore introduce concise notation for the masses of 
the charm mesons and the pions. We denote the masses of the spin-0 charm mesons D° 
and by Mq and Mi, respectively. We denote the masses of the spin-1 charm mesons 
D*^ and D*^ by Af*o and Af*i, respectively. We denote the masses of the pions 7r° and 
7r+ by mo and mi, respectively. (The numerical subscript is the absolute value of the 
electric charge of the meson.) The pion mass scale corresponding to either mo or mi will 
be denoted by m^. The result of a recent precision measurement of the mass by the 
CLEO Collaboration is Mq = 1864.85 ± 0.18 MeV, where we have combined the errors in 
quadrature [11]. We use the PDG values for the other masses [10]. The errors on the pion 
masses are negligible compared to those on the charm meson masses. Some of the differences 
between the charm meson masses have errors that are significantly smaller than the errors 
in the masses themselves. 

We also introduce concise notations for simple combinations of the masses. We denote 
the reduced mass of a spin-1 charm meson and a spin-0 charm meson by 

We denote the reduced mass of a pion and a spin-0 charm meson by 

"^'^ - m, + Mi 

We denote the differences between the D* masses and Dn thresholds by 

5ijk = M^i - Mj - ruk. (4) 

The differences between the D* masses and the thresholds for Dn states with the same 
electric charge are 

5ooo = 7.14 ±0.07 MeV, (5a) 
Son = -2.23 ±0.12 MeV, (5b) 
Sioi = 5.85 ± 0.01 MeV, (5c) 
duo = 5.66 ± 0.10 MeV. (5d) 
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The isospin splittings between the charm meson masses are Mi — Mq ^ 4.8 MeV and 
M,i - M,o ~ 3.3 MeV. The energy sphtting ly = (M*i + Mi) - (M*o + Mq) between the 
D*+D- and D*°D^ thresholds is 

i/ = 8.08 ±0.12 MeV. (6) 

A phenomenological analysis of the decays of the D*^ and D*"*" was presented in Ref. [26] . 
We summarize here the results of that analysis, which was based on chiral symmetry and 
isospin symmetry. The PDG value for the total width of the D*+ is r[£)*+] = 96 ± 22 keV 
[10]. Using the PDG values for the branching fractions for D*~^ decays, we obtain measured 
values for the partial widths for decays: 



T[D*+ D°7r+] = 65.0 ± 14.9 keV, (7a) 

T[D*+ D+7r°] = 29.5 ± 6.8 keV, (7b) 

r[D*+ D+-f ] = 1.5 ± 0.5 keV. (7c) 

Using isospin symmetry and the PDG values for the branching fractions for D*^ decays, we 
obtain predictions for the partial widths for D*^ decays: 

Y[D*° D°7r°] = 40.5 ± 9.3 keV, (8a) 

r[D*o ^ D%] = 25.0 ± 6.2 keV. (8b) 



The prediction for the total width of the D*" is T[D*^] = 65.5 ± 15.4 keV. 

The decay rates for D* Dtt are fairly sensitive to the mass of the D*, since they scale 
like the 3/2 power of the energy difference between the D* mass and the Dn threshold. A 
virtual D*'^ (or £)*+) with energy M*o + E (or M*i + E) can be considered as a, D* whose 
rest energy differs from its physical mass by the energy E. The width of the virtual particle 
varies with E. We denote the energy- dependent widths of the D*^ and D*^ by r^:i{E) and 
T^o{E), respectively. If \E\ is small compared to myr, these energy-dependent widths can be 
obtained simply by scaling the physical partial widths for the decays D* Dn: 

r,o{E) = r[i?*° ^ D%] + r[D*' ^ D%'] [ [{Sooo + E)/6ooof' 0{6ooo + E) 

+2 (mii/moo)'/' [(<5oii + E)l8o^of^ e{5on + E)] , (9a) 

r,i{E) = r[D*+ ^ L>+7] + r[D*+ ^ L>+7r°] [(5iio + E)/5i,of' 0{5no + E) 

+r[L'*+ ^ L>°7r+] [(5ioi + £;)/(5ioi]'/' ^(5ioi + E). (9b) 

We ignore any energy dependence of the decay widths into D7, because the photon energy 
and the phase space for the decays D* — > do not vary significantly in the D*D threshold 
region. In Fig. 1, we plot the energy- dependent widths V^q{E) and r*i(-E' — u) as functions 
of E. The offset z/ 8.1 MeV in T^i{E — u) was chosen so that r^,o(-E') and r*i(£' — u) are 
the relevant widths for a. D*D system consisting of D and a D* with total energy E relative 
to the threshold. Thus r,o(^) reduces to T[D*^] at E ^ and r^i{E - ly) reduces 

to r[D*+] at E ^ u. The physical widths r[D*'^] and r[D*+] are shown in Fig. 1 as data 
points with error bars. At the 0**^0^ threshold, the energy-dependent width of the D*'^ is 
r,i(-z/) ^ 1.5 MeV. 

The individual terms in Eqs. (9) have obvious interpretations as energy-dependent partial 
widths for decays of D*+ and We can define energy-dependent branching fractions by 
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FIG. 1: The energy-dependent widths T^o{E) and r*i(£;— i^) for a virtual D*^ with energy M^o + E 

and a virtual D*~^ with energy M*i + E — v, respectively, as functions of E. The points with error 
bars at £" = and E = u indicate the central values and uncertainties of the physical widths of 
D*^ and D*~^ , respectively. 



dividing these terms by r*i(£') or V^q{E). For example, the energy-dependent branching 
fractions for D*° D^n^ and D*+ D+t:^ are 

Brooo(^) = ^^^'^ 7^?"^"^ 1(^^000 + E)/5,^of^ ^(^ooo + E), (10a) 

i *oi-C'J 

Brno(^) = T^?^"""^ [(^110 + E)/6^,,f' 9{6uo + E). (10b) 

The standard isospin multiplets for the charm mesons are (—D'^,D'^), (D^,D~), 
(— D*^), and (Z)*°, D*^), where the first and second states are the upper and lower com- 
ponents of the multiplet, respectively. The D*D channels with charge conjugation quantum 
number C = + are 

{D*D)l = +j^{D*^D" + D^D*^) , (11a) 
{D*D)l = {D*+D- + D+D*-) . (lib) 

The superscript i on {D*D)\ is the absolute value of the electric charge of either meson. We 
will refer to [D*D)^j^ and (D*Z))^ as the neutral and charged charm meson channels, respec- 
tively. The channels with isospin quantTim numbers / = and 1=1 are the antisymmetric 
and symmetric linear combinations of these neutral and charged channels, respectively: 

= T2 [(^*^)+ - ' (12a) 

iD*DY=' = ^ [{D*D)l + {D*D)\] . (12b) 
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III. LOW-ENERGY D*D SCATTERING 



In this section, we discuss the low-energy scattering of charm mesons. We first summarize 
the results of Ref. [26], which took into account only the neutral channel (£)*L>)° . These 
results should be accurate when the energy E is within a few MeV of the D*^D^ threshold. 
We then extend the region of validity to the entire D*D threshold region by taking into 
account the resonant coupling to the charged channel 



A. Neutral channel only 

We begin by writing down a general expression for the transition amplitude for S-wave 
scattering in the {D*D)^_^ channel that is compatible with unitary. The transition ampli- 
tude A{E) for the scattering of nonrelativistically normalized charm mesons in the channel 
can be written in the form 

A{E) = -^f{E), (13) 

where f{E) is the conventional nonrelativistic scattering amplitude expressed as a function 
of the total energy of the charm mesons. An expression for the scattering amplitude that is 
compatible with unitarity is 

where k{E) — {—2M^qqE — ieyl'^ and E is the total energy relative to the D*^D° threshold 
in the center-of-mass frame. If the inverse scattering length 7 is complex, the imaginary 
part of the scattering amplitude in Eq. (14) is 

Imf(E)^\f{E)\Hm[^-K{E)]. (15) 

The scattering amplitude f{E) in Eq. (14) satisfies the constraints of unitarity for a 
single-channel system exactly provided 7 is a real function of E. For positive real values of 
the energy E, Eq. (15) is simply the optical theorem for this single-channel system: 



lmf{E) = \fiE)\W2M,ooE (E > 0). (16) 

The left side is the imaginary part of the T-matrix element for elastic scattering in the 
[D*D)^ channel multiplied by M*oo/(2vr). The right side is the cross section for elastic 
scattering multiplied by (2M^oo-£')^^^/(47r). We first consider the case 7 > 0. In this case, 
the amphtude f{E) has a pole at a negative value of the energy E, indicating the existence 
of a stable bound state. If 7 varies sufficiently slowly with E that it can be approximated 
by a constant, the pole is near -Epoic ~ — 7^/(2M*oo) and the binding energy is 7^/(2M*oo)- 
In addition to the contribution to the imaginary part of f{E) in Eq. (16), there is a delta- 
function contribution at = ^^poie^ 

Imf{E) ^ ^S{E + jy{2M,oo)) {E < 0,7 > 0). (17) 

-'"*00 

We next consider the case 7 < 0. In this case, the pole in the amplitude f{E) is not on the 
real E axis, but on the second sheet of the complex variable E. The standard terminology 
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for such a pole is a virtual state. The imaginary part of the amphtude is nonzero only in 
the positive E region and is given by Eq. (16). 

Scattering in the (D*D)*]_ channel cannot be exactly unitary, because the D*^ has a 
nonzero width and because the charm mesons have inelastic scattering channels. The in- 
elastic channels include D^D^tt^ and D^D^'-f, which are related to D*^ or D*^ decays, as 
well as all other decay modes of X(3872), including J/-?/' tt+tt" , J/ip7i~^7i~7i'^, and J/ipj- 
In Ref. [26], the dominant effects of the D*^ width and the inelastic scattering channels 
were taken into account through simple modifications of the variables 7 and k{E) in the 
scattering amplitude f{E) in Eq. (14). The effects of the decays of the constituent D*^ 
or were taken into account simply by replacing the mass M*o that is implicit in the 
energy E measured from the D*^D^ threshold by M*o — ir*o(-E')/2, where T^,q{E) is the 
energy-dependent width of the D*^ given in Eq. (9a). This changes the energy variable 
k{E) = (-2M*oo^ - ief''^ into 

k{E) = V-2M*oo[£^ + ir,o(£^)/2]. (18) 

At the threshold E — the energy-dependent width r*o(-E') reduces to the physical width 
r[D*°]. The expression for k{E) in Eq. (18) requires a choice of branch cut for the square 
root. If E is real, an explicit expression for k,{E) that corresponds to the appropriate choice 
of branch cut can be obtained by using the identity 



^J-2M[E + iT/2\ = \fM 



(19) 



In Ref. [26], the effects of inelastic scattering channels for the charm mesons other than 
D^D^Ti^ and D^D^'-f were taken into account by replacing the real parameter 7 by a complex 
parameter with a positive imaginary part. The expression for the imaginary part of the 
amplitude f{E) in Eq. (15) can now be interpreted as the optical theorem for a multi-channel 
system consisting of {D*D)'\_ and all the inelastic scattering channels. The right side can 
be interpreted as the total cross section for scattering in the (D*D)^ channel multiplied 
by (2M*oo-E)^^^/(47r). The terms proportional to ImK(E) and Im7 are proportional to the 
elastic and inelastic cross sections, respectively. This interpretation requires Im7 > 0. 

The scattering amplitude f{E) in Eq. (14) with k{E) given by Eq. (18) and a complex 
parameter 7 has a pole at an energy E'poie that is not on the real axis. If the difference between 
r*o(-E'poie) and r*o(0) is small compared to r[D*°], the pole energy can be approximated by 

E,,,,^--^-zr[D*']/2. (20) 
The energy £^poie is that of a bound state if Re7 > and that of a virtual state if Re7 < 0. 



B. Coupled neutral and charged channels 

We now generahze the results of Section III A to the system consisting of the two coupled 

channels (D*Z))1_ and {D*D)\_ defined by Eqs. (11). The amplitudes for transitions between 
these channels can be expressed as a 2 x 2 matrix Aij{E), i,j G {0,1}. We first write 
down a general expression for the transition amplitudes for S-wave scattering in the two 
channels that is compatible with unitarity in this two-channel system. A convenient way to 
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parametrize these amplitudes is to express the inverse of the matrix of amphtudes Aij{E) 
in the form 

A(F)-^^ — (''^^ ° \(-1O0 + k(E) -701 \ f ^/^Uo \ 

^ ' 27r V Vm:^i J \ -701 -7ii + «i(^)yV VM:^! J ' 

(21) 

where k{E) = {-2M^qqE - iefl'^ , ki{E) = {-2M^u{E — v) — ie)^/^, and E is the energy 
relative to the D*^D^ threshold. The parametrization of the inverse matrix in Eq. (21) 
was chosen so that the analytic expressions for the entries of Aij{E) would be as simple as 
possible. It is convenient to define a matrix fij{E) of scattering amplitudes by 

A,iE) = U,{E). (22) 

The entries of the matrix fij{E) are 

foo{E) = f-7oo + ^{E) X^Ti^l , (23a) 

f„,^E) ^ + 1-7.0 + >cW][-Tn+>..W] r' _ 

V 7oi / 

ME) = (-^n + K,{E) ^^^riWs) ■ ^23c) 

V -7oo + /«(-£/) y 

If the parameters 700, 7oi, and 711 are complex, the imaginary parts of the scattering 
amplitudes in Eq. (23) satisfy 

Im/,,(E) = UE)f;,iE)lm[^oo-f^m + UE)f;,{E)lm[^u-n,{E)] 

+ {UE) f;,{E) + fa{E) f;o{E))lm [701] • (24) 

Since Im fij{E) is real, an alternative form for this unitarity equation can be obtained by 
taking the complex conjugate of the right side. 

The amplitudes fij{E) in Eqs. (23) satisfy the constraints of unitarity for this two-channel 
system exactly if 700, 701, and 711 are all real functions of E. For positive real values of the 
energy E, the expressions in Eqs. (24) for the imaginary parts of the amplitudes foo{E) and 
fn{E) are just the optical theorems for this two-channel system: 

Im/oo(£;) = \foo{E)\W^M,ooE + \foi{E)\W2M,u{E - i^) e{E - v), (25a) 
Im/n(E) = |/oi(£;)|V2M,ooS+|/n(S)|V2M,n(^-/^)^^(£-/^). (25b) 

The left sides of Eqs. (25a) and (25b) are proportional to the imaginary parts of the T-matrix 
elements for elastic scattering in the [D*D)^^ and [D*D)\^ channels, respectively. The first 
and second terms on the right side of each equation are proportional to the cross sections 
for scattering into the (Z}*l))*j_ and [D*DY_^ channels, respectively. In the region £" < 0, the 
imaginary parts of foo{E) and fii{E) may also have delta function contributions analogous 
to the one in Eq. (17). 

There are two hmits in which the amplitudes foi{E) and fu{E) go to and foo{E) 
reduces to the single-channel amplitude in Eq. (14). The first hmit is — > -|-oo, which 
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corresponds to increasing the energy gap between the two thresholds. In this case, foo{E) 
reduces to Eq. (14) with 7 = 700- The second hmit is 701, 711 — > 00 with 701/711 fixed, which 
corresponds to decreasing the interaction strength between the two channels. In this case, 
/oo(-E) again reduces to Eq. (14) but with 7 = 700 — 7oi/7ii- 

Scattering in the {D*D)'^ and {D*D)\_ channels cannot be exactly unitary, because the 
D*^ and D*^ have nonzero widths and because the charm mesons have inelastic scattering 
channels. The inelastic channels include DDtt and DD'j, which are related to D* or D* 
decays, as well as all the decay modes of X(3872), which include J/'^ tt+tt", J/-^ 7r"'"7r~7r°, 
and J/V'7- The effects of decays of D*^ and D*^ can be taken into account by replacing 
k(E) in the amplitudes in Eqs. (23) by the expression in Eq. (18). Similarly the effects of 
decays of D*+ and D*~ can be taken into account by replacing ki{E) in the amplitudes in 
Eqs. (23) by 

Ki{E) = ^/-2M,n[E -v + iT,i{E - iy)/2]. (26) 

At the D*~^D'^ threshold E = i>, the energy-dependent width T^:i[E — i>) reduces to the 
physical width F If E is real, an exphcit expression for k,i{E) that corresponds to the 
appropriate choice of the square-root branch cut in Eq. (26) can be obtained by using the 
identity in Eq. (19). The effects of inelastic scattering channels other than DDtt and DD'j 
can be taken into account by replacing the real parameters 700, 7oi, and 711 by complex 
parameters with positive imaginary parts. The expression in Eq. (24) for the imaginary 
part of the amplitude fii{E) can now be interpreted as the optical theorem for the multi- 
channel system consisting of {D*D)^, {D*D)^, and all the inelastic scattering channels. The 
right side can be interpreted as the total cross section for scattering in the {D*D)\_ channel 
multiplied by (2M*jj£^)^/^/(47r). The terms proportional to ImK(E) and ImKi(E) are the 
cross sections for scattering into the {D*D)^ and {D*D)\_ channels. The terms proportional 
to Im7oo, Im7oi, and Im7ii give the remaining inelastic cross sections. This interpretation 
requires Im7oo > and Im7ii > 0. 

C. Constraints from isospin symmetry 

We now proceed to exploit the approximate isospin symmetry of QCD. Deviations from 
isospin symmetry can be treated as small perturbations except at low energies that are 
comparable to the isospin splittings between hadron masses, which in the case of charm 
hadrons are less than 5 MeV. In strong interaction processes, isospin-symmetry-violating 
effects come primarily from hadron mass differences. Exact isospin symmetry would require 
the masses of the charged charm mesons to be equal to those of their neutral counterparts, 
which implies i/ — and M^n = M*oo- It would also require the inverse matrix of amplitudes 
in Eq. (21) to be diagonal in the isospin basis. These conditions can be expressed as 

rr ( -700 + «(^) -701 \tt\ - ( -7o + «(^) \ , . 

^ \ -701 -711 + <E) )^ -\ -71 + ^{E) ) ' ^^^> 

where 70 and 71 are the inverse scattering lengths in the isospin-symmetry limit for the 
7 = and I — 1 channels, respectively, and U is the unitary matrix associated with the 
transformation between the charged/neutral basis in Eqs. (11) and the isospin basis in 
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Eqs. (12). The conditions in Eq. (27) imply 

700 = (71 + 7o)/2, (28a) 

701 = (7i-7o)/2, (28b) 
711 = (7i+7o)/2. (28c) 

The constraints on the amphtudes fij{E) from the approximate isospin symmetry of QCD 
are obtained by inserting these values for the parameters into Eqs. (23). 

In terms of the parameters 70 and 71, the scattering amplitudes in Eqs. (23) reduce to 

ME) ^ -'^""^;^^'-'^> , (290) 

where the denominator is 

D{E) = 27170 - (71 + 70) [i^{E) + Ki{E)] + 2ki{E)k{E). (30) 

The unitarity conditions in Eq. (24) can be written 

lmf,j{E) = U{E)f;,{E)lm[^^+^o-2n{E)]/2 

+U(E) f;,(E) Im [71 + 70 - 2k,{E)] /2 

+ {fio{E) f*i{E) + fa{E) f;o{E))lm [71 - 70] /2. (31) 

If there is a bound state or virtual state near the D*^D^ threshold with complex energy 
E'poie, the denominator D{E) given in Eq. (30) vanishes at that energy. If we define a variable 
7 by 

7 = K(^pole), (32) 

the equation D{Epoie) = can be expressed as 

-fKi{Epole) - i(7l + 70) [7 + «l(^pole)] + 7i7o = 0. (33) 

The variable 7 can be identified with the inverse scattering length introduced in Eq. (14). 
The energy i?poie is given approximately by Eq. (20). If we neglect the small difference 
between Ki(£'poie) and Ki(0), one can obtain an approximate solution of Eq. (33) for 70 in 
terms of 71 and 7: 

^ 7i^i(0) + 7i7-2/.i(0)7 
' 27i-Ki(0)-7 ■ ^ ^ 

If the energy E is within a few MeV of the D*^D^ threshold, the scattering amplitudes 
in Eqs. (29) can be simplified. If the small diffrence between ki{E) and ki(0) is neglected, 
the denominator D{E) given in Eq. (30) reduces to 

D{E) - [71 + 70 - 2«:i(0)] [-7 + k{E)] . (35) 
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In the numerators, k,{E) and 7 can be neglected compared to ki(0), 70, and 71. The 
scattering amphtudes then reduce to 



foo{E) «i f{E), (36a) 

ho{E) ^ f{E), (36b) 

71 + 7o - 2k,i(0) 

^ ^ /'^ol (36c) 

7i + 7o - 2/ti(0) 

where /(-E) is the single-channel scattering amplitude in Eq. (14). If 7 is neglected compared 
to fi:i(0) and 71, the expression for 70 in Eq. (34) reduces to 

71^1 (0) .oyN 

Using this expression to eliminate 70 in favor of 71, the coefficient of f{E) in the amplitudes 
fij{E) in Eqs. (36) can be factored into a term that depends on the channel i and a term 
that depends on the channel j: 

f,,{E)^cd-{E)cj, (38) 
where the coefficients Cj are given by 

Co — I, (39a) 



The values of the two independent parameters 70 and 71 could be calculated using poten- 
tial models for heavy mesons with pion-exchange interactions. As pointed out by Tornqvist 
in 1993, these models indicate that there should be D*D bound states near threshold in 
several / = channels, including the S-wave l^"*" channel, but not in any of the 1 = 1 
channels [27]. Tornqvist could not predict whether the I'*'"'' state was just barely bound 
or not quite bound, because his results depended on an ultraviolet cutoff whose value was 
estimated to be the same as the corresponding ultraviolet cutoff for the two-nucleon system 
[27]. He also could not predict whether the state would be closer to the D*^D^ threshold 
or the D*^D~ threshold, because his calculations were carried out in the isospin symmetry 
limit. With the discovery of the X(3872), the ambiguity associated with the ultraviolet 
cutoff can be removed by using the observed binding energy of the X(3872) to tune the 
value of the ultraviolet cutoff. One can then use the meson potential model to predict the 
binding energies of other heavy meson molecules in both the charm sector and the bottom 
sector [23]. 

In the absence of explicit calculations of the parameters 70 and 71, one can still use results 
of the meson potential model calculations in Ref. [27] to get some idea of the likely values of 
these parameters. The bound state near threshold with / = and J^'^ = 1+"*" arises from the 
effects of coupled S-wave and D-wave channels. In the S-wave channel, the pion-exchange 
potential is not deep enough to give a bound state. The D-wave interaction provides just 
enough additional attraction to obtain a bound state very near threshold. Thus we expect 
I70I to be significantly smaller than the natural scale m^. The sign of 70 could be either 
positive or negative. The meson potential model calculations in Ref. [27] indicate that there 
IS no bound state with 7 = 1 and J^^ = 1++. The pion-exchange potential has the opposite 
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sign as in the 1 = case, so it is repulsive. We therefore expect 71 to be positive and 
comparable to or larger than the natural scale m.^. In particular, 71 should be much larger 
than |7o|. Given an estimate of 71, an estimate of 70 is actually superfluous because it can 
be determined using Eq. (34). 

The scattering amplitudes fij{E) in Eqs. (29) simplify if the parameter 71 is assumed to 
be large compared to Ki(0). The denominator D{E) given in Eq. (30) reduces to 

D{E)^ -^,[-2^0 + k{E) + ^,{E)]. (40) 

The scattering amplitudes reduce to 



The matrix projects onto the 1 = channel. The denominator in Eq. (40) vanishes at -Epoie- 
If the small difference between Ki(i?poie) and ki{0) is neglected, we get an approximate 
expression for 70 in terms of the variable 7 defined by Eq. (32): 

7.«^. (42) 

In Ref. [25], the authors analyzed data from the Belle and Babar Collaborations on the 
energy distributions of J/t/jn'^n' and D^D^n^ near the X(3872) resonance produced by the 
decay — > + X. Their model for the {D*D)^ elastic scattering amplitude f{E) is a 
generalization of the Flatte parametrization for a near-threshold resonance [28] : 

/hkkn(^) = _(2/^) [E-Ef + ir{E)/2] + ki{E) + k{E) ' ^^^^ 

where k{E) = {-2M^qoE - ief^ and ki{E) = {-2M^oq{E - v) - ief^ . The function T{E) 
is determined up to normahzation factors by the decays of X(3872). The other adjustable 
parameters are g and Ef. In Ref. [25], the authors found that their fits had a scahng behavior 
that made it impossible to determine unique values of the parameters. As pointed out in 
Ref. [22], the scaling behavior simply indicates that their fits were insensitive to the term 
— {2/g)E in the denominator in Eq. (43). If this term is deleted and if —{2/g)Ef and T[E)/ g 
are identified with the real and imaginary parts of 270, the scattering amplitude in Eq. (43) 
reduces to the amplitude foo{E) in Eq. (41), except that it does not take into account the 
effects of the D* widths. As pointed out in Ref. [26], the D*^ width can be taken into 
account by replacing k,{E) by the expression in Eq. (18). Similarly, the D*~^ width can be 
taken into account by replacing Ki{E) by the expression in Eq. (26). 



IV. LINE SHAPES OF X(3872) 

If a set of particles C has total quantum numbers that are compatible with those of 
the X(3872) resonance and if the total energy E of these particles can be near the D*^D^ 
threshold, then there can be a resonant enhancement in the channel C. The line shape of 
X(3872) in the channel C is the differential rate for producing the particles C as a function of 
their total energy E. In Ref. [20], it was pointed out that the line shapes of the X(3872) can 
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be factored into short- distance factors that arc insensitive to E and the inverse scattering 
length 7 and a long-distance factor that is determined by E and 7. In Ref. [22], it was shown 
that the factorization formulas could be derived using the operator product expansion for an 
effective field theory that describes the cc sector of QCD near the D*^D^ threshold. There 
is a factorization associated with the creation of the charm mesons if all the particles in the 
initial state and if the particles in the final state other than the resonating particles in C have 
momenta in the resonance rest frame that are of order 171^^ or larger. If C is a short- distance 
decay mode of X, there is also a factorization associated with the inelastic scattering of 
the charm mesons into the particles in C. A short- distance decay mode of X(3872) is one 
for which all the particles have momenta that are of order or larger in the resonance 
rest frame. Examples of short-distance decay modes are J/ip 7r~^7r~ and J/iJjTT~^7r~7r^. An 
example of a decay mode that is not short-distance is D^D'^n'^. 

In this section, we consider the fine shapes in the decays B K+C, where C is a channel 
that is enhanced by the X(3872) resonance. We first summarize the results of Ref. [26] in 
which only the neutral channel {D*D)^ defined in Eq. (11a) was taken into account. These 
results should be accurate when the energy E is within a few MeV of the D*^D^ threshold. 
We then extend the region of validity to the entire D*D threshold region by taking into 
account the resonant couphng to the charged channel {D*D)\_ defined in Eq. (lib). 



A. Neutral channel only 

Expressions for the line shapes of the X(3872) that take into account the D*^ width and 
inelastic scattering in the {D*D)'^ channel were derived in Ref. [26]. We give here a more 
exphcit derivation of the line shapes produced by the decay S+ — > K'^ + X. Our starting 
point is the optical theorem for the width of the 

r[i3+] = -^lmA[B+ ^ (44) 

where ^[i?"*" — > B~^] is the one-meson-irreducible forward amplitude for B^. This ampli- 
tude has contributions from intermediate states consisting of a K'^ recoiling against sets 
of particles whose invariant mass A/^o + ^0 + E is near the D*^D^ threshold. There is 
resonant enhancement for small E if the particles are accessible from the (D*Z))1_ channel. 
The resonant contributions to the forward amplitude can be expressed as a loop integral 
over the 4-momentum of the K'^: 



^^[s. - «i - y [cr. HE) ci; j (45) 

There is an implicit restriction of the integral to the region of small E. The expression 
inside the parentheses takes into account the amplitude for the creation of charm mesons 
in the channel (D*i))^, the resonant propagation of the pair of charm mesons, and the 
amplitude for their annihilation. Factorization has been used to express it as the product of 
a long-distance factor and two short-distance factors. The long-distance factor f{E) is the 
scattering amplitude for elastic scattering in the {D*D)^^ channel given in Eq. (14). The 
short- distance factors depend on the 4- momenta Pb and Pk of the B^ and , but 
they are insensitive to the small energy E defined by (Pb — PrY = {M^q + Mq + E^. The 
short-distance factors can therefore be simphfied by setting E — 0. Using the Cutkosky 
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cutting rules, the resonant contribution to the imaginary part of the forward amphtude can 
be written 

ImAes[i^+ ^B^]--J {^bI IrnfiE) (Cf^)*) . (46) 

Again there is an imphed restriction of the integral to the region of small E. The contribution 

to the width of 5+ from its decay into and the X(3872) resonance can be obtained by 
inserting Eq. (46) into Eq. (44). The distribution in the invariant mass M — M^q + Mq + E 
of the resonance can be obtained by inserting the identity 

1 = J d'Pn 5\Pb -Pk- Pr) I dM' 5{M' - P^). (47) 

Changing the order of integration and using Eb — Ek > and \E\ <^ M*o + M), this can 
be written 

(48) 

Upon inserting this into Eq. (46) , we obtain a factorization formula for the inclusive energy 
distribution summed over all resonant channels: 

— [S+ -^K+ + resonant] = 2 Tf + Im / (E) . (49) 
dE 



The short-distance factor is a positive real constant: 

M*o + Mo f d^Pk f d^Px /o N4r4/p p jD \\rK+\2 



(50) 



A more explicit expression for the short-distance factor can be obtained by using Lorentz 
invariance to express the short-distance factor in the form 



Cb+ — Pb ■ (cd*)*, (51) 

where e^j* is a polarization vector for the D*^ [18, 19] and is a constant with dimensions 
of inverse mass. Evaluating the phase space integral and summing over the D*^ spins, we 
get 



A3/2(MB,mK,M,o + Mo) 



B+ — 



647r2(M,o + Mo)M| 



2 



(52) 



The optical theorem in Eq. (15) can be used to resolve the inclusive resonant rate in 
Eq. (49) into two terms according to whether they have Im7 or ImK(E) as a factor. We 
interpret the term proportional to Im7 as the contribution from all short-distance decay 
channels C. The imaginary part of 7 can be expressed as a sum over those decay channels: 

Im7 = J]r^(£;). (53) 

c 

We have allowed for the possibility that the dependence of some of the short-distance factors 
r*" (£^) on the energy E may not be negligible in the D*^D^ threshold region. Thus the energy 
distribution in a specific short-distance channel C can be expressed as 

^K^ + C] = 2rS \f{E)\'r^{E). (54) 
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We interpret the term in Eq. (49) proportional to lmK{E) as the contribution from channels 
that correspond to D*^D^ or D^D*^ followed by the decay of the D*^ or Z)*°. We can 
resolve this term into the contributions from the channels D^D'^tt^, D~^D^7r~ , D^D~7r~^, and 
D°L'°7 by multiplying it by the energy-dependent branching fractions Broool-E"), |Broii(£'), 
|Broii(£'), and Broo7(£'), which add up to 1. A simple expression for lmK{E) can be obtained 
by using the identity in Eq. (19). The resulting expression for the energy distribution in the 
D^D^Ti^ channel is 

^[5+ ^ K+ + D'D\'] = 2rf: \f{E)\' [m^ooWE' + r,o{Ey/4 + E)] ^'^ Brooo(i^), 

(55) 

where Brooo(-E') is given in Eq. (10a). 

The energy distributions of the X resonance in the decays + X are given by 

expressions identical to those in Eqs. (49), (54), and (55) except that the short- distance 
constant is replaced by F^o- Thus the line shapes for X(3872) produced in decays 
and B^ decays are predicted to be identical in the region within a few MeV of the D*°D^ 
threshold. 

In Ref. [21], the decay rates of X into J/ip plus tt+tt^, tt'^ti'tt^, n^'-f, and 7 were calculated 
under the assumption that these decays proceed through couplings of the X to J/ip and 
the vector mesons and u. The results of Ref. [21] were used in Ref. [26] to calculate 
the dependence of the factor r'^{E) in Eq. (54) on the energy £■ for C = J/t/jn'^n' and 
J / ip 7r~^ 7r~ 71^ . The normalization factors r'^(O) can only be determined by measurements 
of X(3872) decays. Simple analytic approximations to {E)/r^{0) that are accurate to 
within 1% in the region \E\ < 8.5 MeV for J/ipir^ir" and in the region \E\ < 1 MeV for 
J/V^TT+TT-TT^ are given in Ref. [26]. 

In Figs. 2 and 3, we illustrate the line shapes for X(3872) near the D*^D^ threshold. We 
take into account the D*^ width, but we neglect the effect on the line shapes of inelastic 
scattering channels for the charm mesons. We show the line shapes for three values of 7: 
-F34, 0, and -34 MeV. For 7 = -^34 MeV, the peak of the resonance is at £; = -0.6 MeV, 
which is the central value of the measurement in Eq. (1). In Fig. 2, we show the line shapes 
in a short-distance decay mode, such as J / ip 7r~^ tt^ or J / ip tt~ . The line shape is given 
by Eq. (54). We have neglected the energy-dependence of the factor V'~"{E). The relative 
normalizations of the curves for the three values of 7 are determined by using the same short- 
distance factors F^+ and F'^. For 7 = -|-34 MeV, which corresponds to a bound state, the 
line shape is dominated by the Breit-Wigner resonance near E — —0.6 MeV. For 7 = —34 
MeV, which corresponds to a virtual state, the line shape has a cusp near E = {) MeV. In 
Fig. 3, we show the line shapes in the D'^D^tt^ channel. The line shape is given by Eq. (55). 
The relative normalizations of the curves for the three values of 7 are determined by using 
the same short-distance factor F^+ . For 7 = -|-34 MeV, which corresponds to a bound state, 
the dominant features of the line shape are a Breit-Wigner resonance near E = —0.6 MeV 
and a threshold enhancement for E > 0. For 7 = —34 MeV, which corresponds to a virtual 
state, the hne shape has only the threshold enhancement. 



B. Coupled neutral and charged channels 

We proceed to generalize the factorization formulas in Section IV A to the two-channel 
case. We begin by generalizing the forward amplitude in Eq. (45). We have to take into 
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E (MeV) 

FIG. 2: The line shapes near the D*'^D^ threshold for X(3872) decaying into a short-distance 
channel, such as J/ip-K'^TT^ or J / il) it~^ t:~ tt^ . The line shapes are shown for three values of 7: +34 
MeV (solid line), (dotted line), and -34 MeV (dashed line). 



account the possibility of resonant scattering between any pair of the charged and neutral 
channels. The amplitude can be written 

- 51 - - / g E E (CS:- ME) cr.-) pr^^- (56) 

Following the same path as in Section IV A, we ultimately arrive at a factorization formula 
for the inclusive energy distribution summed over all resonant channels: 

_ [5+ ^ x+ + resonant] = 2 J] J] vf^^^^ Im {E) . (57) 
The short-distance factors are 



The short-distance factors Fg^'^'^ and F^+'^^ are positive real constants, while F^ 



(F^+'^°)* is a complex constant. Thus there are four independent real constants associated 
with the transitions. These constants satisfy the Schwarz inequality 



pJs:+,oi 

^ B+ 



'<rj:'°°F5:'-. (59) 
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I-H 






E (MeV) 

FIG. 3: The line shapes near the 0*^3^^ threshold for X(3872) in the D^D'^tt'^ channel. The line 
shapes are shown for three values of 7: +34 MeV (solid line), (dotted line), and —34 MeV (dashed 
line) . 



The optical theorem in Eq. (31) can be used to resolve the inclusive resonant rate in 
Eq. (57) into four terms according to whether they have Im7o, Im7i, Imn^E), or lmK,i{E) 
as a factor. We interpret the terms proportional to Im7o and Im7i as the contributions 
from short-distance decay channels C. The imaginary parts of 70 and 71 can be expressed 
as sums over those decay channels: 



im7o = J]r^(£;), 

c 

im7i = J]rf(£;). 



(60a) 
(60b) 



The factorization formula for the energy distribution in a specific short-distance decay 
channel C is 



1 1 



E E ^m'nUE) - UE)][f;,{E) - f;,{E)] T^{E) 



1=0 j=0 

1 1 



+ E E ^nt'^iUE) + UmUUE) + fUE)] rf (£;). (61) 

\i=0 j=0 / 

The terms in Eq. (57) proportional to Im/t(i?) and lmK,i{E) also have simple interpreta- 
tions. We interpret the term proportional to lmK{E) as the contribution from channels that 
correspond to D*^D^ or D^D*^ followed by the decay of the D*^ or D*^. We can resolve this 
term into the contributions from the individual channels D^D^n^, D+D^tt", D^D'n'^, and 
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D^D^^ by multiplying it by the energy- dependent branching fractions Broool-E"); |Broii(i?), 
and |Broii(-E), and Broo7(-E), which add up to 1. For example, the line shape of X in the 
D^D^Ti^ decay mode is 



V i=o j=o ' 

V , -1 1/2 

X [m,oo(V^2 + (r.o(^)/2)2 + E)J Brooo(^), (62) 



where Brooo(-E') is given in Eq. (10a). We interpret the term in Eq. (57) proportional to 

\m.Ki{E) as the contribution from channels that correspond to D*~^D~ or D~^D*~ followed 
by the decay of the D*^ or D*~ . We can resolve this term into the contributions from 
the individual channels D~^D~7r^, D^D~7r^, D+i)°7r~, and D^D^'-f by multiplying it by 
energy-dependent branching fractions. For example, the line shape of X in the D'^D~7r° 
decay mode is 



dE 



\ i—n n—n / 



i=0 j=0 

, -1 1/2 

M,n W{E - vY + (r,i(£; - i/)/2)2 + E - v) Brno(£^), (63) 



where Briio(-E) is given in Eq. (10b). The expressions for the line shapes in the decay 
channels D^D~tt^ and D~D^'k^ are more complicated because they receive contributions 
from channels that correspond to D*^D° or D^D*^ as well as channels that correspond to 
D*+D- or D+D*-. 

If the energy E is very close to the D*^D^ threshold, the two-channel factorization formu- 
las in Eqs. (57), (61), and (62) should reduce to the single- channel factorization formulas in 
Eqs. (49), (54), and (55). For the factorization formula for S+ K+ + D'^D'^tt'^ in Eq. (62), 
this can be verified by inserting the expressions in Eq. (38) for the scattering amplitudes 
fij{E) at small E. The factorization formula reduces to Eq. (55) with the short-distance 
factor given by 

rf:-EErf:''^c.c*. (64) 

Similarly, the factorization formula for — > K'^ -|- C in Eq. (61) reduces to Eq. (54) with 
given by Eq. (64) and r*"(£^) given by 

r^(^) - ^^r^(£^) + ^^rf (^). (65) 

To see that the two-channel factorization formula in Eq. (57) for the inclusive resonant rate 
reduces to Eq. (49), we express the imaginary part of fij{E) in a form that is compatible 
with the Cutkosky cutting rules: 

lmfij{E) = a f{E) (Imc,) + q {l^f{E)) c* + (Imq) f^E) c*. (66) 

Since cq = 1, it has no imaginary part. The expression for ci in Eq. (39b) is a function of 
7i and «;i(0) only. The imaginary part of ki(0) is suppressed relative to its real part by a 
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factor of r[D*'^]/i'. We expect 71 to have a real part that is comparable to or larger than 
mjr, so the imaginary part of 71 should also be small relative to its real part. Thus the only 
term on the right side of Eq. (66) that is not suppressed is the one with the factor lmf{E). 
Inserting that term into the factorization formula in Eq. (57), we find that it reduces to 
Eq. (49) with rf+ given by Eq. (64). 

The factorization formulas for the energy distributions simplify if the parameter 71 is 
assumed to be large compared to Ki{0). The scattering amplitudes fij{E) in Eq. (29) reduce 
to the expressions in Eq. (41). The two-channel factorization formula in Eqs. (57), (61), and 
(62) all reduce to the single-channel factorization formulas in Eqs. (49), (54), and Eq. (55) 
with the scattering amplitude f{E) replaced by the expression for foo{E) given in Eq. (41). 
By using Eq. (42) to eliminate 70 in favor of 7, the scattering amplitude reduces to 

^^^^ " -^ + KiE)+\{E)-KM- ^^^^ 
The short-distance factor for 5+ — > K'^ transitions reduces to 

rS^EEl-ir^^rS:''- (68) 

The sums project the {D*Dy_^ channels onto isospin 0. The short-distance factor for the 
short-distance decay channel reduces to 

r^(^) ^ 2T^{E). (69) 

The coefficient of T^{E) goes to zero in this limit. Thus the decay of X into final states C 
with total isospin quantum number I — 1, such as J/-^ tt+tt" , are suppressed in the large-71 
limit. 



C. Constraints from isospin symmetry 

We have not yet fully exploited the approximate isospin symmetry of QCD. Since the 
short- distance factors only involve momenta of order ttItt and larger, isospin-violating effects 
can be neglected in these factors. Thus isospin symmetry can be used to constrain the 
short- distance factors. At the quark level, the transitions B — > K + D*D and B ^ K + 
DD* proceed through two operators in the effective weak Hamiltonian: the charged current 
operator 67^^(1 — 75)007^(1 — 75)5 and the neutral current operator 67'*(1 — 75)5 07^(1 — 75)0. 
These operators are both isospin singlets. Thus isospin symmetry is respected by these 
transitions. It can therefore be used to relate the short-distance coefficients '* for the 

B~^ transition to the short-distance coefficients C^o '* for the B^ transition. 

Since B^ and B^ form an isospin doublet and and form an isospin doublet, the 
coefficients C^o'* and C^+'* are related by Clebsch-Gordan coefficients: 

Cgo ' — ~^B+ ' (70a) 
C'' - -C^t''- (70b) 
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This implies that the short-distance constants F^o in the factorization formulas for S° — > 

transitions are related to the corresponding constants in the factorization formulas 

for transitions by 

= Tf,''\ (71a) 

rS?^ = {T^^'r, (71b) 
r^'^^ = rj:'"^ (71c) 

Thus the short-distance constants associated with the — > and — > transitions 
are determined by four independent real constants. 

Isospin symmetry also constrains the short-distance factors V^^E) associated with decays 
of X into short-distance decay modes. It imphes that for a decay channel C with definite 
isospin quantum number 7 = or 7 = 1, only the term with the factor Ti{E) contributes. An 
example of a decay channel with isospin quantum number 7 = is J/ilJ7r~^7r~7r^, assuming 
that the 7r"'"7r~7r'' comes from the decay of a virtual u. An example of a decay channel with 
isospin quantum number 7 = 1 is J / ip 7t~^ 7i~ , assuming that the tt+tt" comes from the decay 
of a virtual p°. We will give the factorization formulas for short-distance decay channels 
with definite isospin quantum number 7 = and 7 = 1 for both B~^ transitions and 

B^ transitions. For a short-distance decay channel C with definite isospin quantum 

number 7 = 0, such as J/ipTr^T:~TT^, the energy distribution in Eq. (61) reduces to 



dE 



+F5:'-|7.-.(^)r)^, (72a) 



|D(E)|2' 

^[73° ^ 7^° + C] = 4 (fS:'°°|7i - ^{E)\^ - 2Re[FS:'°^(7i - ^{E)){^, - ^,{E)) 



-rS:^-K-..(i^)r)^. (72b) 



Dm- 

For a short-distance decay channel C with definite isospin quantum number 7 = 1, such as 
J/ipTT'^Tr'' , the energy distribution in Eq. (61) reduces to 



dr 

dE 



[B+ ^ 7r+ + C] = 4 (Fg;'°°|7o - K,{E)\' + 2Re[F5;'°^(7o - k,{E)){jo - K{E)r] 

+F-:^-|7o-.(^)r)gg|„ (73a) 
^[B'^K' + C] ^ 4 (rj:'°°|7o - k{E)\' + 2Re[F5:'°^(7o - k{E)){^o - K,{E)r] 

+F-:--|7o-..(^)r)^. (73b) 

In Eqs. (72b) and (73b), we have used the isospin symmetry relations in Eqs. (71) to express 

the short- distance coefficients F^o'*"* in terms of F^+'*-'. 

The effects of the charged charm meson channel {D*D)}^ on the line shapes of X(3872) 
in the decays B ^ K + J/t/jn'^n" and B ^ K + J/'07r"'"7r~7r'^ have been discussed recently 
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by Voloshin [29]. Voloshin made conceptual errors by ignoring resonant scattering between 
the {D*D)\ and {D*D)]_ channels and ignoring the constraints of isospin symmetry on the 
transitions B ^ K. In Voloshin's paper, our parameters 70 and 71 are denoted by kq and 
Ki and the analogs of our functions n{E) and Ki{E) are denoted by —ikn and Kc- Voloshin 
took into account the constraints of isospin symmetry associated with the J / ijj 7r"^7r~ and 
J/ip 7i~^7T~'K^ in the final state. His results for the energy distributions can be expressed in 
the form 



dr 

dE 



[B ^ K + J/i/jn+n-Tr'^] 



7i - Ki{E) 



dT 
dE 



DiE) 
7o - Ki{E) 



D{E) 



(74a) 
(74b) 



where k{E) = {-2M^qqE - ief/'^ and Ki{E) = {-2M^u{E - u) - ief/'^. The normalizing 
factors ^•^/''i'f' and ^■^I'i'^ can presumably be different for B^ decays and B^ decays, although 
this was not stated explicitly in Ref. [29]. The line shapes however were predicted to be 
the same for S+ decays and B" decays. Voloshin's results in Eqs. (74) correspond to 
specific choices for the short-distance factors r^+'*-' in our general factorization formulas 
in Eqs. (72) and (73). In the case of B^ decays, his results in Eqs. (74) are consistent 
with our factorization formulas in Eqs. (72a) and (73a) if r^+''^° is the only nonzero short- 
distance factor for the B ^ K transition. In the case of B^ decays, his results in Eqs. (74) 

are consistent with our factorization formulas in Eqs. (72b) and (73b) if rg+'^^ = Fgo'"" 
is the only such nonzero factor. However these conditions for B^ decays and B^ decays 
are inconsistent. Thus Voloshin's results are incompatible with the constraints of isospin 
symmetry associated with the B ^ K transitions. The primary conceptual error in Ref. [29] 
was the assumption that there is a resonance in the amplitude only if the B ^ K transition 
creates the charm mesons in the neutral channel (D*£))° . However there is also a resonant 
contribution coming from the B ^ K transition creating charm mesons in the charged 
channel {D*D)\_ followed by the resonant scattering of the charm mesons into the neutral 
channel. A second conceptual error in Ref. [29] was the failure to take into account the 
constraints of isospin symmetry on the amplitudes for the B ^ K transition. 



D. Current-current factorization and heavy-quark symmetry 

In Ref. [19], it was pointed out that the combination of a standard current- current fac- 
torization approximation and heavy-quark symmetry could be used to simplify the factor- 
ization formulas associated with the X resonance va. B ^ K transitions. In the standard 
current-current factorization approximation, the matrix elements of the relevant terms in 
the effective weak Hamiltonian are expressed as products of matrix elements of currents: 

{KD*D\hY{l - 75)cc7^(l - i,)s\B) ^ {D\h^{l - ^,)c\B) {KD*\cl^.{l - l,)s\Q), 

(75a) 

{KD*D\h>'{l - 75)5 C7^(l - l^)c\B) ^ {K\h>'{l - iMB) (^*5|c7m(1 - 75)c|0). 

(75b) 

The D* and D in the final state can equally well be replaced by D and D* . The matrix 
element of the charged current 67^(1—75)0 in Eq. (75a) is nonzero only if the D contains the 
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same light quark as the B. In the case of a B^, the D or D* must be or D*^. In the case 
of a -B°, the Z) or D* must be D~ or D*^ . As pointed out in Ref. [19], heavy-quark symmetry 
implies that the matrix element of the neutral current C7^(l — 75)0 in Eq. (75b) vanishes at 
the D*D threshold. Thus this matrix element is suppressed in the D*D threshold region. 
Putting these two observations together, we conclude that the current-current factorization 
approximation together with heavy quark symmetry puts strong constriants on the matrix 
elements of the effective weak Hamiltonian. It implies that in B~^ — > transitions, the 
formation of the X(SS72) resonance is dominated by the creation of charm mesons at short 
distances in the neutral channel {D*D)^. Similarly, in S° — > K'^ transitions, the formation 
of the X(3872) resonance is dominated by the creation of charm mesons at short distances 
in the charged channel {D*D)\. These statements imply that the short-distance coefficients 

This suppression leads to a 
K transitions in the factorization 



= —C^o arc suppressed relative to = — C^o 

hierarchy in the short-distance factors associated with B 
formulas: 



^^+,01 1 
B+ I 



pX+,00 



(76) 

If we assume that Fg^'^^ and |r^+'°^| are negligible compared to r^+'°°, the expressions 
for the line shapes of X(3&72) become rather simple. For a short- distance decay channel C 
with definite isospin quantum number / = 0, such as J/'07r"'"7r~7r°, the energy distributions 
in Eqs. (72) reduce to 



dE 



K+ + C] 
• X° + C] 



4r 



K+,00 

B+ 



AT 



K+,m 

B+ 



D{E) 
71 - <E) 



D{E) 



For a short-distance decay channel C with definite isospin quantum number / 
J/ipTi~^7i~, the energy distributions in Eqs. (73) reduce to 



dE^ 



K+ + C] ^ 4F 



K+,00 

B+ 



7o - Ki{E) 



4F 



K+,00 

B+ 



D{E) 
7o - k{E) 



rf(£;), 



D{E) 



rf(£;). 



For the D^D^tt^ channel, the energy distribution in Eqs. (62) from the B 
and its analog from the 5° — > transition reduce to 



(77a) 

(77b) 
1, such as 

(78a) 

(78b) 
transition 



dE 



K+,00 

B+ 



7i+7o-2ki(E) 



X 



DiE) 

M,oo{VE^ + r,o{E)y4: + E) 



1/2 



Brooo(£^), (79a) 



dE 



[B' 



2r 



X 



K+,00 

B+ 



71 



70 



D{E) 



nl/2 



M,ooWE^ + T,oiE)yA + E) Brooo(^). (79b) 



Note that the line shapes in Eqs. (77), (78), and (79) arc determined by the parameters 70 
and 7i or, equivalently, 7 and 71. The relative normalizations of the rates from the B^ — > K'^ 
transition and from the S+ — > K'^ transition are also determined by 7 and 71. 
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FIG. 4: The line shapes in the D*D threshold region for X(3872) produced by a B'^ — > or 
transition and decaying into J/^ vr+vr^vr'^. The hne shapes are shown for 71 = ±00 and 
three values of 7: +34 MeV (solid line), (dotted line), and —34 MeV (dashed line). 

In Figs. 4, 5, and 6, we illustrate the line shapes in the D*D threshold region for X(3872) 
produced by S — > X transitions. We take into account the D*^ width, but we neglect the 

effect on the line shapes of inelastic scattering channels for the charm mesons. For simplicity, 
we show only the line shapes for the limiting case 71 ±00. Thus the denominators D{E) 
can be approximated by Eq. (40) and numerator factors such as 71 — ni{E) or 71 — n{E) 
can be approximated by 71. The parameters 7 and 70 are related by the pole equation 
270 — 7 — K\{E^o\e) = 0, where Epoie is given in Eq. (20). If we take 70 to be real, then 7 
has an unphysical negative imaginary part. We therefore take 7 to be real and use the pole 
equation to determine the complex parameter 70: 

70 = ^ (A/2M,ni^ + (M,n/M,oo)7' + iM,ii{r[D*^] - r,,(-u)) + 7) . (80) 

We show the line shapes for three real values of 7: +34, 0, and —34 MeV. The corresponding 
values of 70 have real parts 82 MeV, 63 MeV, and 48 MeV, respectively. Their imaginary 
parts are all approximately 0.00012 MeV, which is completely neghgible. For 7 = +34 MeV, 
the peak of the resonance is at = —0.6 MeV, which is the central value of the measurement 
in Eq. (1). 

In Fig. 4, we show the line shapes in the short- distance decay mode J/ipTT^7T~TT^. The 
line shapes, which are the same for X produced by a — > or B° — > transition, 
are given in Eqs. (77). The relative normalizations of the curves for the three values of 7 
are determined by using the same short-distance factors r^+'°'' and Y'^^^'"'^'" . In Fig. 5, 
we show the line shapes in the short-distance decay mode J/z/'tt^tt". The line shapes are 
given in Eqs. (78). The upper and lower panels show the line shapes produced by — > 
and B° — > K° transitions, respectively. The hne shapes from the S+ — > K'^ transition have 
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FIG. 5: The line shapes in the D*D threshold region for X(3872) produced hy a B ^ K transition 
and decaying into J/tp-K^-K~. The line shapes are different for X produced by a — >■ 
transition (upper panel) and a B'^ ^ K'^ transition (lower panel). The line shapes are shown for 
7i = iboo and three values of 7: +34 MeV (solid lines), (dotted lines), and —34 MeV (dashed 
lines) . 
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E (MeV) 

FIG. 6: The line shapes in the D*D threshold region for X(3872) produced by a or 
transition and decaying into D^D^tt^. The line shapes are shown for 71 = ±00 and three 
values of 7: +34 MeV (solid line), (dotted line), and —34 MeV (dashed line). 

approximate zeros near +6 MeV, while the line shapes from the B^^ — >■ transition have 
approximate zeros near —2 MeV. The relative normalizations of all six curves are determined 
by using the same short-distance factors r^+'''° and Y'^^^'^'^'^ . In Fig. 6, w^c show the line 
shapes in D^D^n^. The line shapes, which are the same for X produced by a or 
B^ transition, are given in Eqs. (79). The relative normalizations of the curves for 

the three values of 7 are determined by using the same short-distance factor V ^j^'^^ . 

Fig. 5 illustrates the fact that the hne shape of the X(3872) may depend not only on 
the decay channel but also on the production mechanism for the resonance. The difference 
between the line shapes in the J/ip-K^n' decay channel for X produced by B^ and 
B^ — > transitions is particularly dramatic because of the approximate zeros in the line 
shapes. These approximate zeros are general features of the line shapes in Eqs. (78). If the 
imaginary parts of 70 and ki{E) are neglected, the numerator factor I70 — ki{E)\'^ in the 
energy distribution in Eq. (78a) has a zero between the D*'^D^ and D*^D~ thresholds. If 
I7I ^ |'ti(0)| ^ |7i|, the approximate expression for 70 in Eq. (34) reduces to /ti(0)/2. The 
zero is therefore near |i/ 6.1 MeV. If the imaginary parts of 70 and k{E) are neglected, 
the numerator factor I70 — /t(-E')P in the energy distribution in Eq. (78b) has a zero below 
the D*^D^ threshold. If I7I < |fi;i(0)| < |7i|, the zero is near --^v ^ -2.0 MeV. In the 
case of B^ decays, the approximate zero forces the line shape to be narrower on the trailing 
edge of the resonance. In the case of B^ decays, the approximate zero forces the line shape 
to be narrower on the leading edge of the resonance. 

In Ref. [19], Braaten and Kusunoki predicted that the decay rate for B^ — > -\-X{^%12) 
should be suppressed relative to that for B^ — > ir+-|-X(3872). Their prediction was based on 
the current-current factorization approximation and heavy-quark symmetry. Together they 
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imply that, in the D*D threshold region, the — > transition creates charm mesons 
predominantly in the neutral channel {D*D)^^, while the — > transition creates them 
predominantly in the charged channel (Z)*D)^. Since the X(3872) is a resonance in the 
channel, the authors of Ref. [19] concluded that the rate for B^ ^ -\-X must be 
suppressed relative to that for B^ — > -\- X. In retrospect, this prediction was the result 
of a conceptual error. 

The conclusion of Ref. [19] that B^ — > + X is suppressed follows from the factorization 
formulas in Sections IV B if ki{Q) is assumed to be much greater than 70 and 71. In this 
hmit, the two-channel factorization formulas for B^ decays in Eqs. (77a), (78a), and (79a) 
all reduce to the single-channel factorization formulas in Eqs. (54) and (55), where f{E) is 
the single-channel scattering amplitude in Eq. (14), 7 = (71 +7o)/2, T{E) = 2Ti{E), and 

— r^+'°°. For the inverse scattering length 7 to be small compared to 70 and 71, 70 and 
7i must be nearly equal in magnitude but opposite in sign. If 7 is small compared to 71, the 
factorization formula for S° decays in Eqs. (77b), (78b), and (79b) also reduce to the single- 
channel factorization formulas in Eqs. (54) and (55), with = [7i/Ki(0)]^r^+ . The rate 
for B^ + X is therefore suppressed by a factor of [7i/ki(0)]^ compared to the rate for 

B^ — > + X. Thus the conclusion of Ref. [19] is consistent with the factorization formulas 
only if the parameters satisfy the hierarchy I7I <^ |7i| <^ |/ti(0)|. Since the conclusion of 
Ref. [19] is only valid in one corner of the parameter space, the authors must have made a 
conceptual error. 

The conceptual error in Ref. [19] has to do with the momentum scale at which the in- 
ferences from the current-current factorization approximation and heavy-quark symmetry 
are applied. Heavy quark symmetry is relevant at energy scales that are small compared 
to the heavy quark mass rric and large compared to the energy scale set by isospin symme- 
try violations, which is v ^ 8.1 MeV. Equivalcntly, it is relevant at momentum scales that 
are small compared to ttIc and large compared to ki(0) ~ 125 MeV. The current-current 
factorization approximation and heavy quark symmetry imply that the short-distance con- 
stant Fg^'™ = r^o'^^ dominates. This inference should be applied at a momentum scale 
where heavy quark symmetry applies, which requires the momentum to be large compared 
to fi:i(0) ~ 125 MeV. The conceptual error in Ref. [19] was inferring that the transitions 
B^ — > + [D*D)^_^ and B^ — > + {D*D)^ dominate at a momentum scale small com- 
pared to «;i(0). At this low momentum scale, there is a resonance only in the {D*D)^_^_ 

channel. However the dominance of F^+ '*^° = F^o'^^ at a scale large compared to ki(0) does 
not imply the dominance of the transitions 5+ -\- {D*D)\ and B^ ^ + {D*D)\ at 

lower scales. As the momentum scale is lowered, resonant scattering between the {D*D)'^ 
and {D*D)l channels can feed the transitions B+ K+ + {D*D)\ and 5° ^ K^ + {D*D)\. 

To deduce the correct implications of the current-current factorization approximation and 
heavy quark symmetry at momentum scales small compared to «;i(0), we can consider the 
general factorization formulas in Eqs. (77), (78), and (79) in the low energy region where 
/t(-E) is small compared to Ki(0). In this region, the factorization formulas reduce to the 
single-channel factorization formulas in Eqs. (54) and (55). The assumption that F^+'°° = 
F^o'^^ dominates implies that the short-distance constants for the B ^ K transition are 
Fg+ fti F^+'°° and F^o k, F^+'°° |cip, where the coefficient ci is given in Eq. (39b). The 
ratio of the rates is |cip = |7iP/|7i — /«i(0)p. If we ignore the small imaginary parts of 71 
and ft:i(0), this ratio is greater than 1 if 71 > «;i(0)/2 and less than 1 if 71 < ki(0)/2. Thus 
the rate for B^ ^ -\- X need not be suppressed compared to that for B^ -\- X. 
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V. SUMMARY 



In Ref. [26], we derived line shapes of the X(3872) that should be accurate in the region 
within a few MeV of the D*^D^ threshold. The hne shapes were derived from an expression 
for the resonant scattering amplitude in the {D*D)^ channel that takes into account the D*^ 
width and inelastic charm meson scattering channels. In the factorization formulas for the 
line shapes, short- distance effects and long-distance effects are separated into multiplicative 
factors. The line shapes of Ref. [26] are independent of the production mechanism for the 
X resonance. The line shape in D^D^n^ is different from the line shape in a short-distance 
decay mode, such as J/V'tt+tt^tt^ or J/ip7i~^7i~. As shown by the analysis of Ref. [26], 
the difference in these line shape can explain the difference between the masses of X(3872) 
measured in the J/ipTT~^'!T~ and D^D^ti^ decay modes [6, 7]. 

In this paper, we have derived hne shapes for the X(3872) whose region of validity extends 
to the entire D*D threshold region by taking into account the resonant coupling between the 
(D*Z))'j_ and (D*Z))^ channels. By taking into account isospin symmetry at high energies, 
the coupled-channel scattering amplitudes were expressed in terms of two parameters: the 
7 = and I — 1 inverse scattering amplitudes 70 and 71. Isospin symmetry was also 
taken into account in the short-distance factors in the factorization formulas. In the case 
of production of the X resonance m B ^ K transitions, isospin symmetry reduces the 
short-distance factors to four independent real constants: r^^'^*^, r^+'^^, and the real and 

imaginary parts of r^+'°^. The resulting factorization formulas for the inclusive resonance 
production rate in the — ^ transition is given in Eq. (57). The factorization formula 
for the D^D^n^ channel is given in Eq. (62). The factorization formulas for 7 = and 7 = 1 
short-distance decay channels arc given for both the B^ and the B^ — >■ transitions 

in Eqs. (72) and (73). The line shape in an 7 = short-distance decay channel, such as 
J/^/) 7r"'"7r~7r°, is different from the line shape in an 7 = 1 short-distance decay channel, such 
as JZ-^TT+TT". The line shapes for the X resonance produced by the 5+ — > transition 
are also different from the line shapes produced by the B^ — > transition. 

If we use the current-current factorization approximation together with heavy quark 
symmetry, the factorization formulas simplify dramatically. The factorization formulas for 
an 7 = decay channel, an 7 = 1 decay channel, and D^D^ir^ are given in Eqs. (77), 
(78), and (79), respectively. The short-distance constants associated with the B ^ K 
transitions reduce to a single real constant r^+'°'^. Thus the ratios of production rate 
in 5+ — > K'^ transitions and in B^ — > transitions are completely determined by the 
scattering parameters 70 and 71. 

Our results allow us to identify conceptual errors in previous work on this problem. In 
Ref. [19], Braaten and Kusunoki predicted that B^ ^ -\- X should be suppressed by at 
least an order of magnitude compared to B^ -\- X. The prediction was based on 

the current-current approximation and heavy quark symmetry. The conceptual error was 
an implicit assumption that 70 and 71 are small compared to «;i(0). In Ref. [29], Voloshin 
predicted that the hnes shapes for the X resonance produced by B^ — > -\- X and 
B^ ^ -\- X should be identical. The conceptual errors were ignoring resonant scattering 
between the (7}*_D)° and {D*D)\ channels and failing to take into account isospin symmetry 
in the B ^ K transitions. 

Our results provide a physical interpretation for the model of the (7^*7))° scattering 
amphtude used in Ref. [25] . The scaling behavior of the fits in Ref. [25] indicate that the 
term —{2/g)E in the inverse of the scattering amplitude in Eq. (43) can be omitted. The 
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resulting scattering amplitude is essentially equivalent to the general (D*D)'^ scattering 
amplitude foo{E) in Eq. (29a) in the hmit |7i| ^ |/ti(0)|. 

The most important parameters for predicting the line shapes are the scattering param- 
eters 7o and 71. They could be determined phenomenologically from ratios of rates for 
+ X and + X. Alternatively they could be calculated using the meson 

potential model of Ref. [27]. If these scattering parameters were calculated, the predictive 
power of the results of this paper would be dramatically increased. The line shapes for 
the X(3872) resonance produced hy B ^ K transitions also depend on the short-distance 
factors r^+'°°, r^+'°^, and r^+'^^. They could be determined phenomenologically from 
measurements of the charm meson invariant mass distributions in the decays B — > K + D*D 
and B ^ K + DD*. 

The accuracy of our predictions for the line shapes could be further improved by taking 
into account pions exphcitly. The system consisting of D*^D^, D^D*^, and D^D^n^ states 
with energies near the D^D*^ threshold can be described by a nonrelativistic effective field 
theory. The simplest such theory has S-wave scattering in the {D*D)'^ channel and tt*^ 
couplings that allow the decay D*^ — * D^n^. Fleming, Kusunoki, Mehen, and van Kolck 
developed power-counting rules for this effective field theory and showed that the pion 
couphngs can be treated perturbatively [30] . They used the effective field theory to calculate 
the decay rate for X(3872) D^D^tt^ to next-to-leading order in the pion coupling. In 
applying this effective field theory to the line shapes of the X(3872), one complication that 
will be encountered is infrared singularities at the D*^D^ threshold that are related to the 
decay D*^ — > D°7r°. This problem has been analyzed in a simpler model with spin-0 particles 
and momentum-independent interactions [31]. The problem was solved by a resummation 
of perturbation theory that takes into account the perturbative shift of the D^D*^ threshold 
into the complex energy plane because of the nonzero width of the D*^ . 

In summary, the establishment of the quantum numbers of the X(3872) as and the 
measurement of its mass imply that it is either a charm meson molecule or a charm meson 
virtual state. These two possibilities can be distinguished in practice by their different 
predictions for the line shapes of the X(3872). The analysis of Ref. [26] indicates that 
the existing data favor a charm meson molecule, but a virtual state is not excluded. The 
expressions for the line shapes used in that analysis should be accurate only within a few 
MeV of the D*^D° threshold. The expressions for the hne shapes derived in this paper 
should be accurate in the entire D*D threshold region. When more extensive data on the 
line shapes of the X(3872) in various decay channels and for various production processes 
becomes available, it should be possible to determine conclusively whether the X(3872) is a 
bound state or a virtual state of charm mesons. 
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